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Abstract 

We derive the finite temperature conductance peak distributions and peak-to- 
peak correlations for quantum dots in the Coulomb blockade regime assuming 
the validity of random matrix theory. The distributions are universal, depend- 
ing only on the symmetry class and the temperature measured in units of the 
mean level spacing, A. When the temperature is comparable to A several 
resonances contribute to the same conductance peak and we find significant 
deviations from the previously known T < A distributions. In contrast to 
the T <C A case, these distributions show a strong signature of the charging 

energy and charge quantization on the dot. 
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Quantum dots are two-dimensional microstructures of micron scale or smaller in which 
a small number of electrons are confined by electrostatic potentials. They can be fabricated 
with relatively little intrinsic disorder, in which case the motion of the electrons is ballistic. 
The transport properties of such dots are measured by coupling them to leads through point 
contacts. As the contacts are pinched off, the dot becomes more closed (i.e. more weakly 
coupled to the leads) and the electron resonances become well-isolated. For temperatures 
that are low compared with the mean level spacing, the dot's conductance is dominated 
by the resonance that is closest to the Fermi energy of the electrons in the leads. Since a 
tunneling event requires the addition of one electron into the dot and a collective charging 
energy of e 2 /C (where C is the capacitance of the dot), the conductance exhibits a series 
of approximately equally spaced peaks as a function of the gate voltage 0. The height of 
the conductance peaks, however, shows order of magnitude fluctuations. These fluctuations 
measure directly the fluctuations of the wavefunctions in the interface region between the 
dot and the leads. Using random matrix theory (RMT), the statistical distribution of 
conductance peak heights was derived in closed form 0. Recently, these distributions were 
measured in dots with single- channel symmetric leads and several hundred electrons, for 
both the case of conserved and broken time- reversal symmetry and were found to be 
in good agreement with the theoretical predictions. The measured parametric correlator of 
the conductance peak as a function of an applied magnetic field was also found to be in 
agreement with the predicted correlator 0. 

One aspect of the data from ref. |4| has remained unexplained. Strong correlations 
were observed between the heights of adjacent peaks, in contrast to the RMT prediction of 
vanishing correlations in the low temperature limit. In this experiment the temperature was 
estimated to be 0.3 — 0.5A, so the correlations could be due to deviations from the T < A 
results; but then the rather good agreement of the conductance distributions with the T < A 
result is somewhat puzzling. We address these questions here by deriving conductance peak 
distributions and peak-to-peak correlations for temperatures that are not much smaller than 
the mean-level spacing A. We find that, due to effects of the charging energy, the corrections 
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to the distributions are smaller than expected from the non-interacting Landauer-Biittiker 
conductance formula |J. Nonetheless significant peak-to-peak correlations are induced, and 
we will discuss their relation to experiment below. The deviation of the finite temperature 
distributions from those predicted by naive application of the Landauer-Biittiker formula is 
somewhat surprising since on resonance the mean interaction energy difference between the 
N and N — l particle ground states vanishes J7|J§ and indeed in the limit T<A one obtains 
exactly the same distribution as in the non- interacting case PJlQl- Our result is derived by 
combining the theory of sequential resonant tunneling in the Coulomb blockade regime f7||| 
with the statistical assumptions of RMT. 

Beenakker considered the linear response of a dot in equilibrium with a chemical 
potential equal to the Fermi energy in the leads. The electrostatic energy of the dot with iV 
electrons is given by U(N) = (Ne) 2 /2C — NeaV g where C is the capacitance of the dot, V g 
is the gate voltage and a denotes the ratio between the plunger gate to dot capacitance and 
the total capacitance. The Fermi energy of a resonant tunneling event of the iV-th electron 
into the dot is given by the condition E F = U(N) — U(N — 1) + E N , where E N is the 
single-particle energy of the iV-th level. Alternatively we can define an effective Fermi energy 
E F = E F + eaV g for which the condition of resonant tunneling is E F = E N + (N — 1/2) e 2 /C. 
If the gate voltage is tuned to satisfy eaV g = (N — 1/2) e 2 jC one reaches the degeneracy 
condition for which the charging energy U(N) — U(N — 1) vanishes (to order A). Assuming 
this condition, and that T is always greater than the resonance width, Beenakker showed 
that the resonant conductance G(T, Ep) can be written as a weighted sum over the T = 
resonances A in the dot 

G(T, E F ) = e i§pg where g = Ex w x (T, E F )g x . (1) 

Here g x = 2f _1 [r / A r^/(r / A + T r x )] are the level conductances, where r A is the width of a 
resonance level A to decay into the left (right) lead and T = T l + T r is the total average 
width. The quantity g\ is dimensionless and temperature- independent. w\ = w\(T,Ep) is 
the weight with which a given resonance A contributes to the conductance. The contribution 
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to w\ from any fixed number of electrons on the dot is the product of the probability that 
the level E x is filled with the dot having that number of electrons, and the probability that 
there is an empty state in the leads at the corresponding total energy [|7| : 

•,2 
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P/v is the probability that the dot has N electrons, (n\) N is the canonical occupation of a 
level A, and /(e) = [1 + exp^/fcT)] -1 . In many experiments T, A e 2 /C and only one 
term in @ contributes to a given conductance peak, corresponding to Nq electrons in the 
dot. Eq. (|) reduces to 



w x = 4f(AF No - E F )(n x ) NQ [l - / (E x - E F )\ , (3) 

where AFn = Fjf — Fn-i and Fn is the canonical free energy of N non-interacting particles. 
Here and in the following eaV g is measured relative to (iVo — l/2)e 2 /C 

In the limit T < A only the central level A = N (denoted by A = in the following) 
contributes to a given conductance peak in ([!]). Its weight w becomes the appropriate 
weight for non-interacting electrons, which one would get by appropriate approximation of 
the Landauer-Biittiker formula for narrow isolated resonances 
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In the absence of interactions this result generalizes trivially to the regime where T is 
not much smaller than A. In this case several resonances A contribute to ([I]) with weights 
w^ B obtained by replacing E in (f|) with E x . Since the charging energy "vanishes" on 
resonance one might have expected Eq. @ to reduce to this form. In fact this only happens 
when e 2 /C A; not in the experimentally relevant limit e 2 /C A. If e 2 /C — > 0, then 
all terms (with various number of electrons N) contribute to (0) . The factor 1 — / becomes 
independent of iV and by definition J2n PN(n x ) N = f(E x — Ep) is just the grand-canonical 
occupation number, so that w x reduces to w^ B for all A's. In this case the various manifolds 
of many-electron levels with Nn, No ± 1, AT ± 2; . . . electrons on the dot differ from each 



other only by an energy of order A, and consequently many of the Pjy are non-negligible 
and contribute to W\. However, when the charging energy is large compared with A only 
two manifolds (Nq and Nq — 1) are degenerate while all others are pushed away amounts 
of order e 2 /C. Consequently, the weights w\ differ significantly from their non-interacting 
values when T ~ A. 

We now evaluate these differences quantitatively. This requires the calculation of the 
canonical quantities Fn and (n\) N in (f|), through a projection on a fixed number of particles 
N. This is done in terms of an exact quadrature formula [0 that expresses the canonical 



partition function Zn = e Fn ^ t in terms of grand-canonical partition functions 
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Here the quadrature points are m = 2irm/N sp (N sp is the number of single-particle states), 
and n is a chemical potential chosen anywhere in the range E N < \i < E N+1 . Ei — fi are 
just the particle (i > fi) or hole (i < fi) energies, and cr^ = 1 for a hole and —1 for a particle. 
Since the factors in each term in (|5]) decay exponentially as we move away from /i, only a 
finite number of single-particle states N sp around [i are needed for an exact calculation. The 
canonical occupations are calculated by a similar projection method and differ significantly 
from the corresponding Fermi-Dirac occupations at temperatures of order A or less. They 
lie on a curve similar in shape to the a Fermi-Dirac distribution with a chemical potential of 
/i = (En + E N+ i)/2, but with a level- dependent effective temperature which in the vicinity 
of n is smaller than the actual temperature by almost a factor of 2 jl2| . The inset to Fig. [I] 
shows both distributions for T = 0.5A assuming a uniformly spaced (picket-fence) spectrum. 

While the numerical calculations below include fluctuations of the single-particle energies, 
their effect turns out to be quite small, so that the simple picket-fence spectrum can be used 
to illustrate and understand the results. For a picket-fence spectrum AF N = E N and 
Pn = I(En — Ep). In Fig. [L] we show the weights w\(T,Ep) versus Ep for T = 0.5A for 
several levels around the central level A = 0. The functions w\ become shallower and broader 
as we move away from the level A = 0, in contrast to the shape of the Landauer-Biittiker 



weights (£|), which is independent of A. Assuming the conductance peaks at Ep = E we 
also show in Fig. [i] the weights w\ = w\(T,Eq) that contribute to the conductance peak 
height vs. E x @. 

For the picket-fence spectrum P No = 1/2 and w = (rio) N (denoted in the following 
by (no)); whereas for all levels A ^ 0, the relation w\ ~ w LB /2 holds to within 20% or 
better. Since (n ) < 1 = w LB the actual weights (and hence the conductance) are always 
smaller than predicted by the non-interacting theory. Hence in a rather subtle manner the 
charging energy manifests itself in a suppression of the finite temperature conductance and its 
fluctuations. In the limit T ^> A the weight for the central level Wq = (n ) — > 1/2 = w LB /2 
(see left inset of Fig. []]), and we recover the classical result @ G ~ G LB /2. In fact we find 
that this limit is practically reached at T « 2 A where wq is within 20% of 1/2. However 
a second interesting effect occurs for 0.1 A < T < 2 A. In this interval (no) > 1/2 and the 
ratio w\/w LB is enhanced for the central level relative to adjacent levels. Thus effectively 
the distribution of g is less sensitive to temperature than would be expected from the non- 
interacting theory. 

To test this notion quantitatively we calculated the conductance distributions from Eqs. 
(0)~dl) an d compared them to the non-interacting distributions. A statistical theory of 
the conductance in irregularly shaped quantum dots was developed in Ref. f|| in the limit 
T <C A. The partial width amplitude to decay into a channel c from a resonance level A 
is expressed as the projection of the resonance wavefunction on the channel wavefunction 
across the interface between the dot and the lead. When the electron dynamics in the dot is 
chaotic, the statistical fluctuation of the resonance wavefunction are well described by RMT, 
and the universal distributions of the dimensionless level conductance g\ can be derived. For 



T <C A and single-channel leads, P(g) = y2/irge 2g for conserved time- reversal symmetry 
(GOE) and P(g) = Ag[Ko(2g) + Ki(2g)]e~ 2g for broken time-reversal symmetry (GUE) 




PJI^] (for the case of multi-channel leads see Refs. ]TD|jr5| ). For temperatures comparable 



to A, several resonances contribute to the same conductance peak according to (|I]). The 
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conductance peak is now affected not only by the ffuctuations of the eigenfunctions but also 
by the statistics of the energy levels E\. However we can ignore the fluctuations of the 
energy levels to a good approximation (see below) and take a picket-fence spectrum. The 
conductance peak distribution can then be evaluated in closed form. Assuming the peak is 
positioned at E F = E , the weights w\ are fixed numbers, and only the g x fluctuate. Since 
in RMT different eigenfunctions are uncorrelated we find that the characteristic function of 
the conductance peak distribution P(t) = (e wt ) is given by 



Fig. H shows the analytic distributions derived from (Q) for both conserved (left) and broken 
(right) time-reversal symmetry, and for temperatures T = 0.1A (dotted lines), 0.5A and 
A (solid lines). At the lowest temperature (0.1 A) the distributions essentially coincide 
with those found earlier for T < A, but already at T = 0.5A we see a deviation. The 
non-interacting distributions for T = A are shown for comparison by the dashed lines in 
Fig. |2|; they differ substantially from the actual T = A distributions. Both the mean 
and the variance of the conductance are suppressed at a given temperature due to the 
presence of the charging energy. In the inset to Fig. |2] we compare the measured distribution 
(solid diamonds) for non-zero magnetic field M with the actual distribution we calculate for 
T = 0.3A (solid line), the non-interacting distribution at T = 0.3 A (dashed line) and the 
T <C A distribution (dotted line). The finite temperature distribution describes correctly 
the dip observed in the lowest g data point. The reduced sensitivity to temperature (as 
compared with the non-interacting case) is already observed at this low temperature. 

To test our analytic approximation (|6|), we have done full random matrix simulations 
which include the fluctuations in the single-particle energy levels as well as the possible 
fluctuations in the peak's position. The results are shown by the histograms in Fig. Q The 
largest deviations are observed for T = 0.5A, but even here the analytic approximation 
appears to work well. 




(GOE) 



(GUE) 
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Finally we calculated the peak-to-peak correlator (in a peak sequence vs. gate voltage) 
at finite temperature c(n) = (Gn +tiGn — Gn q 2 ^ / (G 2 No — Gjv f) 2 ), where Gn is the conduc- 
tance peak due to N Q electrons on the dot. In the approximation that the position of each 
peak is fixed at E F = E No + (N — l/2)e 2 /C we can express c{n) in terms of the weights 
W\{Nq) = w\(T,En ) given by (|3|) for No electrons on the dot. Since in both the GOE 
and GUE the eigenvector distribution is independent from the eigenvalues distribution, and 
gxg^ = g^Sxfj, + ~g~x 2 (l — 5a m ) where ~g~\ and g\ are independent of A, we find 



/ x _ Ea w x (N + n)wx(N ) 
c[n) ~ = . (7) 

EA^A(iVo) 2 

The remaining average in (0) is over the energy levels E\, although to a good approximation 
one can again take a picket-fence spectrum. For the latter case we can use the relation 
W\(N + n) = w\-. n (No) to express c(n) w J2x w \-nW\/ J2x w \ m terms of the weights 
w\ = W\(N ) of a fixed number of electrons A^ on the dot. The top inset of Fig. || 
shows c(n) versus n for several temperatures. Fig. ^ itself shows the correlation length 
(defined as the full width at half maximum) as a function of temperature. The increase of 
the correlations between neighboring peaks with T/A is in qualitative agreement with the 
experimental results which show that the peak distributions measured at higher temperature 
|| are more strongly correlated than those measured at lower temperatures 0. The bottom 
inset of Fig. |3] shows a sequence of conductance peaks for a particular realization in RMT at 
T C A and T = 0.5A. Nonetheless, the peak series in Ref. 0] exhibits even stronger peak 
correlations than we would expect for the estimated temperature of T m 0.3 — 0.5A. The 
origin of this enhancement of the correlations at low temperatures is not fully understood 

In conclusion, using RMT we have derived the finite temperature conductance peak 
height distributions and peak-to-peak correlations in Coulomb blockade quantum dots. The 
charging energy was shown to affect the distribution of the conductance peaks at finite tem- 
perature when compared with the results of a non-interacting theory. Further measurements 
of the mesoscopic fluctuations of the conductance at temperatures comparable to the mean 
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level spacing are necessary for a detailed comparison with our theoretical results. 
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FIGURES 




FIG. 1. Left: The weights w\(T,Ep) versus Ep at T = 0.5A and for resonance levels A = 
(solid line), ±1 (dashed lines), ±2 (dotted lines) assuming a picket-fence spectrum. Right: The 
weights w\ (Eq. (|3|), solid squares) at Ep = Eq versus E\ in comparison with the non- interacting 
weights w^ B (solid circles). The right inset shows the canonical occupations (n\) (solid squares) 
versus E\ in comparison with a Fermi-Dirac distribution (dashed line). The left inset shows 

= (^o) as a function of temperature. 
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FIG. 2. The conductance peak height distributions P{g) for the GOE (left) and GUE statistics 
(right). The solid lines are the analytic distributions (see @) at T/A = 0.5,1 for a picket-fence 
spectrum using the weights (|3|), while the histograms describe the respective distributions obtained 
from full RMT simulations (see text). For comparison we also show the T <C A distributions 
(dotted lines) and the non- interacting distributions for T = A (dashed lines). The inset compares 
the experimental distribution in the presence of magnetic field Q (diamonds) with the T = 0.3A 
GUE distribution (solid line), the T = 0.3A non-interacting distribution (dashed) and the T <C A 
distribution (dotted). 
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FIG. 3. The correlation length (full width at half maximum) of the peak-to-peak correlator 
c(n) as a function of temperature. Left inset: c(n) versus peak separation n at several temperatures 
T/A = 0.5, 1, 2. Right inset: a realization of an RMT peak sequence at T <C A (bottom part) and 
T = 0.5A (top part). 
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